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We carry out an analysis of the potential of the Large Hadron Collider (LHC) to discover super- 
symmetry in runs at yfs = 7 TeV with an accumulated luminosity of (0.1-2) fb _1 of data. The 
analysis is done with both minimal supergravity and supergravity) models with nonuniversal soft 
breaking. Benchmarks for early discovery with (0.1-2) fb^ 1 of data are given. We provide an up- 
date of b-tagging efficiencies in PCS 4 appropriate for LHC analyses. A large number of signature 
channels are analyzed and it is shown that each of the models exhibited are discoverable at the 
5cr level or more above the standard model background in several signature channels which would 
provide cross checks for a discovery of supersymmetry. It is shown that some of the benchmarks are 
discoverable with 0.1 fb _1 of data again with detectable signals in several channels. 

PACS numbers: 12.60. Jv,14.80.Ly,95.35.+d 

I. INTRODUCTION 

The Large Hadron Collider (LHC) is currently running and collecting data at yfs = 7 TeV. It is expected that 
it would continue until it has collected at least 1 fb -1 of integrated luminosity. It is then interesting to investigate 
the discovery potential of LHC for supersymmetry (SUSY) in these early runs. 1 Some work in this direction has 
already been done [2HZ] ■ We carry out the analysis of the SUSY discovery potential of the LHC using the next-to- 
lightest supersymmetric partner (NLSP) as a discriminant of models. Several possibilities arise, the most dominant 
of which involve a x i ^? t-i CP odd Higgs A or g as the NLSP. As is well known, the parameter space of minimal 
supersymmetric standard model (MSSM) is very large, consisting of more than a hundred parameters, and thus an 
analysis within this full space of MSSM is intractable. This space is significantly reduced in well motivated models 
such as the minimal supergravity model (mSUGRA) where one assumes R parity invariance and defines the symmetry 
breaking parameters at the grand unification (GUT) scale. These are evolved to low scales by renormalization group 
evolution where the entire sparticle spectrum is determined in terms of the GUT scale parameters. Specifically in 
mSUGRA under the constraint of R parity one has four parameters and one sign [see Eq.(3)] and in nonuniversal 
supergravity models with nonuniversalities in the gaugino sector the space is extended to six parameters and a 
sign. The allowed parameter space is subjected to further theoretical and experimental constraints, i.e., constraints of 
radiative electroweak symmetry breaking (REWSB) with color and charge conservation, satisfaction of flavor changing 
neutral current (FCNC) constraints, and experimental lower limit constraints on the Higgs masses and on the masses 
of the sparticles. Further, in supergravity unified models under the constraints of R parity invariance the lowest R 
parity odd particle (LSP) turns out to be the neutralino over a large part of the parameter space and thus a possible 
candidate for dark matter. Imposition of the WMAP constraint and constraints from direct detection of dark matter 
further reduce the allowed parameter space in these supergravity models. 

However, even after the imposition of all the constraints the allowed parameter space is still large and additional 
criteria must be used in model selection for a study of their signatures. This study is focussed on SUSY discovery 
in the early runs at the LHC, i.e., at y/s — 2 TeV with (1-2) fb" 1 of integrated luminosity. Obviously then it makes 
sense to choose those SUSY models which can be explored within the parameters of the data that will be collected in 
these early runs. However, within this general rubric one ought to make the search as large as possible. For example, 
there is no theoretical necessity to limit our analysis exclusively to just one region of REWSB such as just to the stau 
coannihilation branch. Thus in our analysis, under the constraint that the model points we explore be observable in 
the early runs, we search for a diverse set of models where the NLSP could be one of the allowed set by REWSB. For 
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1 For a review of the discovery potential at y/s = 14 TeV and y/s = 10 TeV see ff] and the references therein. 
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A display of the processes analyzed and their standard model backgrounds at -^/s = 7 TeV 
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TABLE I: An exhibition of the standard model backgrounds computed in this work at y/s = 7 TeV. All processes were 
generated using MadGraph 4.4 |12j . Our notation here is that I = e,/i, r, and all — I, u, jets. Cutsl-Cuts4 indicated in the 
table are defined in Eq. (nj). In the background analysis we eliminate double counting between the process W + tb (ib) and tt 
by subtracting out double resonant diagrams of tt when calculating W ± + tb (ib). 3 



the analysis we present here we find this set to be a chargino (x ± ), a stau (f), a gluino (<?), CP odd Higgs (A ), and 
a stop (t). Thus in the analysis of this paper we give several benchmarks with this diverse set of NSLPs. It should 
be kept in mind that there is still a significant amount of subjectivity in the choice of the benchmarks. To mitigate 
the subjective element in the choice of the benchmarks we further choose the input parameters to be as diverse as 
possible. The investigation is done both for mSUGRA as well as for models with nonuniversalities. The outline of 
the rest of the paper is as follows: In Sec. Q we discuss the analysis of standard model backgrounds at ^fs = 7 
TeV. In Sec. (Ill) we give an analysis of the cross sections for the production of various sparticle processes. These 
include the production of gluinos and squarks, production of a combination of gluinos and charginos or neutralinos, 
and production of charginos and neutralinos. In Sec. (IV) we give an analysis of possible signatures, discuss the SUSY 
discovery potential and give reach plots of the LHC early runs. In Sec. (Jvj) we give several benchmarks c onsi stent 
with all current constraints which are possible candidates for early discovery. Conclusions are given in Sec. (VI) and 
larger tables have been relocated to Sec. (VII). 



II. ANALYSIS OF STANDARD MODEL BACKGROUNDS AT %/s = 7 TEV 

A central element in the discovery of new physics is the determination of the standard model (SM) backgrounds 
for the processes where one expects to see new physics. One such analysis at y/s = 7 TeV has already appeared 
in the literature [!]• Here we give an independent analysis of the relevant backgrounds. In our analysis, we use 
MadGraph 4.4 [T^] for parton level processes, Pythia 6.4 [T3] for hadronization, and PCS 4 for detector simulation [14] . 
We used MLM matching with a kr jet clustering scheme to prevent double counting of final states. Furthermore, 
the 6-tagging efficiency of PGS 4 was updated to better represent the analysis at the LHC. A discussion of this 
improvement is given below. The result of our analysis is presented in Table ^ with parton level cuts as specified in 
Eq. ([T]) and in the caption of Table Our analysis compares well with the analysis of 0]. The differences between 
the two analyses are in part due to the differences between matrix-element Monte Carlo generators which are known to 



2 For some previous works on early discovery though at higher energies see [2 8 11 
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A comparison of b-tagging efficiency in PGS 4 vs ATLAS detectors 

b-tagging Efficiency for E T b-tagging Efficiency for r| 
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FIG. 1: Left panel: A comparison of the fe-tagging efficiency of ATLAS and the loose and tight efficiencies of PGS 4 as a 
function of Et- Right panel: A comparison of the b-tagging efficiency of ATLAS and the loose and tight efficiencies of PGS 4 as 
a function of r\. Ours fits to the efficiency of ATLAS as a function of Et, and r\ as parametrized by Eq. p| are also exhibited. 



exist |15j . For example, we have exclusively used MadGraph with a fc^-based jet clustering algorithm in our analysis, 
while AlpGen [16] uses a cone-based jet clustering algorithm and was the dominant tool used in [4]. In the analysis 
here we have used CTEQ6L1 [T7] parton distribution functions for generating the SM background, and a basic cut 
was applied such that all final state partons (except the top quarks) are required to have > 40 GeV. 

Cutsl = 40 GeV < E T (Jl) < 100 GeV, Cuts2 = 100 GeV < E T (ji) < 200 GeV, Parton level 

Cuts3 = 200 GeV < E T (ji) < 500 GeV, Cuts4 = 500 GeV < E T (Ji) < 3000 GeV cuts. (1) 

An important object in many possible SUSY discovery channels is a b-tagged jet. In PGS 4 b tagging is done based 
off of the Tevatron 6-tagging efficiency. However, it is pertinent to ask if this is a valid approximation of what is to be 
expected at the LHC. Using the Technical Design Reports (TDR) of CMS [H] and of ATLAS T9] one can extract the 
expected b-tagging efficiency as a function of jet pseudorapidity, r\, and transverse energy Et for jets originating from 
heavy- flavor partons. In Fig. we give a comparison of the functions given in the ATLAS TDR [HI [19] with the one 
given in PGS 4. We have omitted the CMS data, since in the TDR b tagging is binned into two sets of rj so one could 
not extract a continuous function from it. The left panel of Fig. ([I]) gives the b-tagging efficiencies as a function of Et 
for ATLAS and the so called "tight" and "loose" efficiencies as defined in PGS 4. One finds a significant difference 
between these and those expected in the ATLAS and CMS detectors. A similar analysis but as a function of r) is 
given in the right panel of Fig. ([!]) where the difference between the efficiencies given by PGS 4 and by (normalized) 
ATLAS is even more glaring. For this analysis, we do not extend the definition of b-tagging beyond r\ of 2.0. Previ 
ously, b-tagging efficiencies were assumed to approach a constant value for Et > 160 GeV. We have extended this 
such that the tagging efficiency reaches a constant for Et > 300 GeV. Thus we have updated 6-tagging functions 
as given in Eq. ^ where we have kept the same degree polynomial as the preexisting ^-tagging functions in PGS 4. 
Furthermore, the total b-tagging efficiency is the product of the Et and r\ functions listed in Eq. Here we make 
no modification to the default PGS 4 rate for mistagging b jets. Our revised b-tagging functions have the form 

b ET = 0.0781391 + 0.0202661£ T - 0.000259664^ + 1.5509 x 10~ 6 £^ - 4.46698 x 10" 9 ^ + 4.7995 x lO" 12 ^ 
b ri = 1.00885- 0.0497485r? + 0.6930367/ 2 - 0.0361142r7 3 - 0.0222204r7 4 + 0.00797621?7 5 . (2) 



III. SPARTICLE PRODUCTION CROSS SECTIONS AT yfs = 7 TEV 

The soft breaking sector of the MSSM is rather large consisting of over a hundred arbitrary parameters. Here we 
use the framework of high scale supergravity to reduce this arbitrariness. Thus, for the supergravity grand unified 



3 When studying W + tb (tb) processes there is a potential to double count such final states if one also considers t, t prodcution processes. 
To prevent this double counting we have eliminated all diagrams involving a top quark from the set of diagrams that lead to W + tb 
final states, with an analogous requirement for W + tb production. 
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Sparticle Production Cross Sections of mSUGRA at ^/s = 7 TeV 
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FIG. 2: An exhibition of the sparticle production cross sections at the LHC at yfs = 7 TeV for mSUGRA as a function of 
the universal gaugino mass mi/j at the GUT scale when mo = 500 GeV, Aq = 0, tan/3 = 20 and sign(^i)= +1. Left panel: 
production cross sections of gg, gq, qq (solid red, dashed green, dashed blue lines). Middle panel: production cross sections 
for gx^, gx° (solid red, dashed green lines). Right panel: production cross sections for X ± X ± i X ± X°) X°X° (solid red, dashed 
green, dashed blue lines). 



model with universal soft breaking, i.e., mSUGRA 20 , one has just four parameters and the sign of the Higgs mixing 
parameter /i; i.e., one has 

mo, mi , Aq, tan/?, sign(/i), (3) 

where m is the universal scalar mass, mi is the universal gaugino mass, A is the coefficient of the trilinear coupling, 
and tan/3 is the ratio of two Higgs vacuum expectation values in the MSSM. Since the nature of physics at the 
Planck scale is still largely unknown, one may extend the minimal supergravity model to include nonuniversalities 
in the gaugino and Higgs sectors as well as in the flavor sector consistent with the FCNC constraints. One of the 
most widely used extensions of mSUGRA consists of the inclusion of nonuniversalities in the gaugino sector [2"lT 
127] . Thus one may include these nonuniversalities by parametrizing the gaugino masses at the grand unified scale, 
which we take to be 10 16 GeV, via the relations m, = m 1 / 2 (l + <5j) (i=l,2,3) corresponding to the gauge groups 
U(l),SU(2) L: andSU{3) C - 

In this framework, the model points are generated by the imposition of REWSB, particle mass limits from LEP 
and the Tevatron, relic density constraints from WMAP [28], the g M — 2 constraints, and FCNC constraints from 
B s — > /i + /i~ and b -> s + 7. WMAP has measured FloMh 2 to a great accuracy with n DM h 2 = 0.1109 ± 0.0056 [2"5] , 
However, to account for the errors in the theoretical computations and possible variations in the computation of the 
relic density using different codes we take a rather wide range in the relic density constraints, i.e., 0.06 < VloAih 2 < 
0.16, in our analysis. Regarding the g^ — 2 constraint [U] recent analysis of the hadronic corrections [30] indicate 
a significant deviation around 3.9<r between the SM prediction and experiment. Such a contribution can arise from 
supersymmetry |31j and the size of the correction indicates the supersymmetric particles (sparticles) to be low in 
mass. 

For the FCNC process B s -> fi + fj,- we take the constraint to be BTZ(B S -> < 5.8 x 10" 8 [31 [33J and for 

the process b — > s-f we take the constraint to be BIZ (b — > sj) = (352 ± 34) x 10 -6 [3U [3S]. We note in passing that 
currently there is a small discrepancy between the SM prediction and the experimental result for b — > 57 which is 
a possible hint for the SUSY contribution |36j and hence another indication of possible low-lying sparticle masses. 
Thus this discrepancy along with the reported — 2 result is encouraging for an early SUSY discovery [37] ■ In 
addition to the above, LEP and Tevatron mass constraints on the sparticle masses and on the Higgs masses are 
applied. These are > 85 GeV, m H ± > 79.3 GeV, m^ > 101.5 GeV, and m fl > 98.8 GeV where A is the 

CP odd Higgs and if is the charged Higgs. Further, we impose the lightest CP even Higgs mass constraint [38] 
m h > (93.5 + 15a; + 54.3x 2 - 48. 4x 3 - 25.7a; 4 + 24.8a; 5 - 0.5) GeV where x = sin 2 (/3 - a) and a is the Higgs mixing 
angle. The final term in the bound represents a theoretical error of 0.5 GeV in the calculation of Mh and Ma assumed 



by the authors of Ref. [35] . Additionally we use the constraints m ± > 104.5 GeV if 



m x ± - m x c 



> 3 GeV for the 



chargino mass and rrig > 309 GeV for the gluino mass [33J. In the analysis we use a top (pole) mass of m t = 173.1 
GeV. Finally, we also apply the constraint that the model points are consistent with recent data from CDMS-II [3j5] and 
XENON-100 [40]. Using MicrOMEGAs 2.4 [41] the spin independent neutralino-proton cross section was calculated 
and compared to CDMS-II and XENON-100. Furthermore, we compare our results to the expected sensitivity 
for XENON-100 6000 kg x day and XENON-lTon for 1 ton x year @2] as well as the expected sensitivity for 
SuperCDMS 03J. 
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Contours of ctsusy in the mg — m x ± mass plane at y/s = 7 TeV 

m = 250 GeV ,tan p= 10 m„ = 250 GeV , tan p = 30 m„ = 1000 GeV , tan p = 10 m = 1000 GeV , tan p = 30 
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FIG. 3: Contour plots with constant values of log(asuSY /fb) for ctsusy in rrig — m x ± mass plane for the case with nonuniversali- 
ties in the gaugino sector. Gaugino masses mi, ma, and 777,3 vary up to 1 TeV. First panel from left: mo = 250 GeV, tan/3 = 10 
while Ao = 0, sign(/i)=+l; second panel from left: mo = 250 GeV, tan/3 = 30; third panel from left: mo = 1000 GeV, 
tan f3 = 10; fourth panel from left: mo = 1000 GeV, tan/3 = 30. 



Next we present the cross sections for sparticle production processes in mSUGRA calculated from Pythia 6.4 at 
y/s = 7 TeV. This is done by generating 5K events for multiple m 1 / 2 values where the other parameters are taken to 
be mo = 500 GeV, Aq = 0, tan j3 = 20, and \x > 0. The results are exhibited in Fig. ^ where the left panel gives the 
cross sections for the production of gg (solid red line), gq (dashed green line), qq (dashed blue line) as a function of 
77ii/ 2 . The middle panel gives the cross sections for the production of gx^ (solid red line), gx° (dashed green line), and 
the right panel gives the production cross section for x ± X ± (solid red line), x ± X° (dashed green line), x°X° (dashed 
blue line). The analysis of Fig. ^ shows these cross sections to be significant, indicating that at low mass scales 
as many as 10 4 or more SUSY events will be generated with 1 fb -1 of integrated luminosity at the LHC. A similar- 
analysis for the case with nonuniversalities in the gaugino mass sector is given in Fig. ([3]), where we give contour plots 
in the m~ g — m x ± mass plane with other parameters as stated in the caption of the figure. The plots give contours of 
constant log {ctsusy /&>) m the range 1 — 3.5. The contours' plots indicate that a chargino mass up to about 500 GeV 
and a gluino mass up to roughly 1 TeV would give up to 10 3 or more events with 1 fb -1 of integrated luminosity. Of 
course, the discovery of sparticles requires an analysis including the backgrounds and selection of appropriate cuts to 
enhance the signal to the background ratio. This will be discussed in Sees. ( IV ) and |v| . 



IV. SIGNATURE ANALYSIS AND SUSY DISCOVERY REACH AT y/s = 7 TEV AND 1 FB 1 

A number of works already exist which analyze the signatures for supersymmetry at y/s — 14 TeV and y/s — 10 TeV 
(for a small sample see [9HTT| [44] ) . In our analysis we computed the sparticle branching ratios with SUSY-HIT [45] 
and they were subsequently read into Pythia. PGS 4 was used for the detector simulation with no trigger imposed 
(L0). Here we give the analysis for y/s = 7 TeV for a large number of signatures for each candidate model considered. 
These are listed in Table |n| and consist of a combination of multijets, b-tagged jets, multileptons, jets and leptons, 
and photons with a variety of cuts geared to reduce the standard model background and enhance the SUSY signal 
with and without missing energy. We discuss first the model with universal soft breaking, i.e., the mSUGRA model. 
In Fig. (J4j) we give the reach of the early runs at the LHC for mSUGRA in the 777,1/2 — m o plane. One finds that the 
reach can extend up to about 400 GeV for 777i/ 2 at low values of m and up to about 2 TeV for m for low values 
of rai/2 with 1 fb -1 of integrated luminosity, and up to about 450 GeV for 2 fb -1 of data, and for mo the reach can 
extend up to 1.9 (2) TeV for l(2)fb _1 of integrated luminosity. For the nonuniversal case a plot in m 1 / 2 — m is not 



very illuminating because of the presence of nonuniversalities. Here we exhibit in Table (III) a set of benchmarks 



which have a chance of early discovery. A criterion used in the selection of these benchmarks is the size of the cross 



section for the production of SUSY events ctsusy which is shown in the last column of Table (III). Here one finds 
that ctsusy for some of the benchmarks is as large as 10-20 pb or more implying that as many as (1 — 2) x 10 4 SUSY 
events will be produced at the LHC with 1 fb -1 of integrated luminosity. Thus with efficient cuts to reduce the 
SM background there appears a good chance for the discovery of such models. A detailed analysis of the signatures 
implementing the cuts of Table |n]) bears this out. Thus as exhibited in Fig. ^ one finds that all of the benchmarks 



of Table (III ) do indeed produce visible signals not just in one but in several channels. In fact for most the benchmark 



models of Table (III) one has as many as five channels and often more where the SUSY signal will become visible, thus 
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Reach Plot at ^/s = 7 TeV up to 2 fb 1 of Integrated Luminosity. 

mSUGRA LHC Reach @ 7 TeV 
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FIG. 4: A reach plot for mSUGRA using the signature analysis and the standard model backgrounds of this work given in 
Table |l]) at the LHC with ^/s = 7 TeV and 1 fb -1 of integrated luminosity. The mSUGRA parameters used are Aq = 0, 
tan/3 = 45, sign(/i) =1. The analysis is done under the conditions of REWSB and the LEP and Tevatron constraints but 
without the imposition of the relic density and FCNC constraints. The condition used for a signal to be observable is S > 
max(5\/ SM, 10) where SM stands for the standard model background. Early LHC reaches at lfb -1 for the gluino (g), the 
chargino {jet), the neutralino (xi)> the stau (fi), and the stop (ti) are exhibited in the inset where the y axis is plotted on a 
logarithmic scale. 



providing important cross-checks for the discovery of supersymmetry. In Fig. (|5| we also exhibit discovery channels 
for 0.1 fb -1 , 1 fb -1 , and 2 fb -1 of data which shows the identities of the new signature channels that open as one 
increases the integrated luminosity. It is interesting to ask how the number of visible signatures depends on the 
integrated luminosity. An analysis of this issue is given in Fig. ^ where a plot is given of the number of signature 
channels where the SUSY signal becomes visible as a function of the integrated luminosity. The figure shows that 
the number of discovery channels increases rather sharply with luminosity and can become as large as 10 or more at 
1 fb -1 of integrated luminosity. Quite interestingly the analysis of Fig. ^ also exhibits that a SUSY discovery can 
occur with an integrated luminosity as low as 0.1 fb -1 still with several available discovery channels. 

As noted earlier all of the benchmarks listed in Table ( III I are discoverable and further, as shown in Fig. ([7| , all 
of them are also consistent with the current limits on the spin independent neutralino-proton cross sections from 
CDMS-II and XENON- 100 and can be seen in the next generation of xenon and germanium experiments. Finally in 
Table (IV) we exhibit the light sparticle spectrum for a subset of the benchmarks given in Table (III). Here we note 



that some of the models in Table (III ), as, for example, the model CI, have typically a small value of /i indicating that 
they reside on or near the so-called hyperbolic branch/focus point region [JS] of radiative breaking of the electroweak 
symmetry. Such models would have scalar masses which are typically heavier and often much heavier than the gaugino 
masses as can be seen in Table ( IV ) . 



V. MODEL SIMULATIONS 



We give now a further discussion of the model simulations. As Fig. ([5| indicates, the primary discovery channels 
for supersymmetry at = 7 TeV will be jet-based signatures which are designed to be as inclusive as possible to 
increase the number of signal events. This will need to take precedence over signal purity (i.e. efficiency to reject 
background) at values of integrated luminosity at or below 1 fb . Four of the five chargino NLSP benchmarks can 
be discovered via jet-based signatures within the first 100 pb -1 of data, with the remainder (Chargino 2) reaching a 
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A Display of Visible Discovery Channels for 0.1 fb -1 , 1 fb _1 and 2 fb" 1 at yfs = 7 TeV. 
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FIG. 5: An exhibition of the visible discovery channels for 0.1 fb 1 (black squares), 1 fb 1 (dark gray squares) and 2 fb 1 
(light gray squares) at y/s = 7 TeV. The discovery channels are listed in Table [IT] Here we are using the convention that C 
denotes a chargino model, G denotes a gluino model, H denotes a Higgs model, S denotes a stau model, and T denotes a stop 
model. 

A Display of the Number of Discovery Channels as a Function of Integrated Luminosity. 



E 



13 
12 
11 
10 



C5 



0.1 



I? 



S6 



L 



G2 



L(fb" 1 ) 



Tl 



FIG. 6: An exhibition of the rapid rise in the number of discovery channels vs integrated luminosity for four early discovery 
benchmarks given in Table (III I. The number of discovery channels for supersymmetry in each case is in excess of five and in 
some cases as large as 10 or above at 1 fb -1 of data at ^fs — 7 TeV. 



five-sigma significance in the multijetslOO channel within 200 pb . The stau NLSP models favor traditional multijet 
signatures such as multijets200 and Ht 500 which involve much harder jet-pT requirements. Among our benchmarks 
the stau NLSP cases tend to have the heaviest gluinos with 709 GeV < rrig < 912 GeV and lightest squarks in the 
range of 330-600 GeV. These examples will therefore exhibit very similar characteristics in the early data taking, with 
the heavier gluino producing long cascades of moderately energetic jets which satisfy the high-pT jet requirements 
common in analyses at s/s = 14 TeV. By contrast, our Higgs, stop, and gluino NLSP models favor discovery channels 
with much looser jet requirements to increase the signal size. These include the Ht 400, multi-bjetsl and multi-bjets2 
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A Display of Benchmarks as Dark Matter Candidates. 
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FIG. 7: An exhibition of the spin independent neutralino-proton cross section, <tsi, for the benchmark models. These are 
labeled by the NLSP, which under the applied constraints allow chargino (green circles), stau (purple triangles), gluino (blue 
squares), CP odd Higgs (red diamond), and stop (orange inverted triangles) NLSPs. In the plot the curve labeled XENON100* 
is the expected sensitivity of XENON100 with 6000kg x days of data, the curve labeled XENON1T is the expected sensitivity 
for lton x year of data and the curves labeled SuperCDMS25 and SuperCDMSIT are the expected sensitivities for the two 
SuperCDMS experiments. 



channels. We can understand the prevalence of b-jet signatures because of the rather small mass gaps between the 
lightest SU (3) -charged state (i.e. the gluino or squark) and the LSP which eventually appears at the end of the 
cascade [57]. For example, the mass gap between the gluino and the LSP for the gluino NLSP models ranges from 60 
to 160 GeV, with an average gluino mass for these models of 390 GeV. We note that these b-jet channels are defined 
with a requirement that pij, > 40 GeV for the b-tagged jet and we impose a veto on isolated leptons (e and /i). The 
latter is not strictly necessary to achieve a high signal significance, but imposing the leptonic veto sufficiently reduced 
the SM background arising from the t,i background samples to increase signal significance by a factor of about 65% 
on average. 

The discussion above indicates an important concern for experimentalists facing the challenge of extracting the most 
signal possible from the first year of LHC data. Traditional signature definitions are often designed to enhance the 
signal-to-background purity so as to be able to make meaningful exclusive measurements of key quantities - often the 
edges or end points of various kinematic distributions. As we will see below, there will likely be few opportunities to 
employ these analysis techniques with only 1 fb _1 of data at ^/s — 7 TeV, even if we are fortunate enough to discover 
supersymmetry in multiple channels. It thus may make sense to seek, at the expense of purity, a more inclusive set 
of signature definitions. Consider, for example, the object defined by 



When the four objects entering Q are restricted to the four hardest jets, the above is typically referred to as effective 
mass (M ff), and has often been studied in the context of multijet channels for discovery of supersymmetry [481 149) . 
When we wish to restrict our attention to jets only, we will refer to Q as the "jet effective mass" of the event. A more 
inclusive definition is to allow the four hardest visible objects to enter Q - in particular, hard leptons. For the sake 
of clarity, when we wish to consider this more expansive variable we will refer to it as "generalized effective mass" or 
Ht, as we did in Table [TT| In the case of this generalized effective mass we will also impose an overall collective cut 
on the scalar sum of all Pt's for these objects, as opposed to cuts on individual object pt's. In an environment where 
jet-pT measurements are likely to be less accurate than we might like, this is a reasonable variable to employ and we 
find it to be particularly effective across many of our benchmark categories. Similarly, we have given in Table pi] ) a 
number of variants to the jet-p-r requirements to give the reader the sense of how relaxing these requirements may 
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Jet Effective Mass and Missing Transverse Energy Distributions 
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FfG. 8: Left: Effective mass distribution (multi-jetslOO) of the four hardest jets for three benchmark models with strong signals 
in this channel: Chargino 1 red (heavy) line, Chargino 5 green (thin) line and Stau 6 blue (dash-dotted) line. This signature 
is defined in Eq. Q in the text. Signal curves are signal + SM background, with the background component indicated by the 
hatched region. Right: Missing transverse energy distribution for the same three models on a logarithmic scale. Again the SM 
component is indicated by the hatched region. 



accelerate discovery. We will investigate this further with a few of our benchmarks below. In general, the lighter jet 
requirements tend to result in a signal which becomes significant sooner, but when both channels are significant the 
stronger jet case may often have the better signal-to-background figure. 

We find very few discovery channels in the leptonic sector for our benchmark models. In the one-lepton plus jets 
channels at I fb _1 wc find between 10 and 50 events for our chargino NLSP and stau NLSP models and less than 10 
events for the Higgs, gluino and stop NLSP models. This is to be contrasted with 367 background events with the 
40 GeV jet requirement (l-lepton40) and 350 with the stricter jet requirement (1-Ieptonl00). These channels achieve 
five-sigma significance only for model Chargino 5, and only after 2 fb _1 of data. The opposite-sign (OS) dilepton 
channels give O(10) signal events after 1 fb _1 of data for our chargino NLSP models only. These cases just fail to reach 
five-sigma significance until 2 fb _1 of data. Event rates for same-sign dileptons and trileptons are consistent with 
zero for all benchmarks at 1 fb _1 , with a handful of events expected at 2 fb _1 . These "gold-plated" supersymmetry 
channels will need to wait for higher center of mass energies or larger data sets to make their presence manifest. 

Finally, we make a few remarks about discovery channels with an imperfect detector. Early results from both 
detectors at the LHC seem to indicate that measurement of missing transverse energy appears to fit Monte Carlo 
predictions within reason and the overall error in determining missing Et is less severe than many pessimists had 
feared [50] - This is welcome news, as our benchmarks will all require missing Et as part of any search strategy. None 
of the benchmarks we studied are visible above backgrounds without this important ingredient, with the number of 
signal events becoming comparable to the SM background only when $ T > 200 GeV or higher is employed. We also 
point out that many of our best signals involve tagged b jets. As mentioned previously, we updated the PGS4 b tagging 
to reflect the more optimistic expectations of the ATLAS and CMS collaborations. Achieving these efficiency targets 
early in LHC data taking will be important for many supersymmetric models relevant to the LHC at y/s = 7 TeV. 

To understand these features we will look a bit more closely at eight of our benchmarks from Table (III) chosen 
to represent an array of phenomenologies at the LHC. These eight benchmarks are collected in Table (IV) where 
we provide ph ysi cal masses of some of the key superpartners relevant for LHC signatures. For each benchmark we 
give in Table |v]) the number of signal events in 1 fb _1 for a selection of signatures from Table |ll|, as well as the 
standard model background count and the signal significance. In this table we have included, for illustrative purposes, 
a modification of the multi-jet signature in which we merely require at least four jets, each with Pj? > 40 GeV in 
addition to the missing transverse energy and transverse sphericity cut. The table illustrates the relative paucity 
of leptonic signatures across all of our benchmark models and the importance of employing relative loose jet-pr 
requirements in multijet signatures. 

In the left panel of Fig. ^ we exhibit the multijetslOO signature, defined by Eq. Q, for the three benchmarks 
of Table |V| with the strongest signal. The distribution of m e g- for the three models is given in terms of the signal 
plus the SM background, normalized to 1 fb _1 , with the SM background also shown separately in the hatched region. 
The strongest signal comes from Chargino 1 whose overall distribution is very similar in size and shape to that of 
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Changing Jet-px Requirements on Multi-Jet Signatures 
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FIG. 9: The effect of changing the minimum transverse momentum of leading jets for multijet signatures is exhibited for 
two models: Chargino 1 (left panel) and Stau 6 (right panel). Heavy (red line) distributions represent the strict multijets200 
signature, thin (green line) distributions are the softer multijetslOO signature, and medium (blue line) distributions are for the 
multijets40 signature. The legend gives the minimum p 1 ^ requirement for the four hardest jets in the event. In this figure we 
show the signal superimposed on the SM background separately to distinguish between the three cases. 



the SM background, while Chargino 5 and Stau 6 peak at slightly larger values. These distributions are heavily 
influenced by the distributions in $ Tl which are plotted on a logarithmic scale in the right panel of Fig. (j8j). We can 
understand these shapes from the values in Table ( IV ) : the relatively light gluino of model Chargino 1 is the lightest 
S'[/(3)-charged state which will decay to low-pr jets, while the lightest such state in the other two cases is the stop 
which has a slightly higher mass (m^ = 402 GeV for Chargino 5 and 497 GeV for Stau 6). For these two cases the 
dominant supersymmetric production processes involve the associated production of a gluino and a squark, with the 
decay of the gluino to on-shell squarks producing slightly harder jets. An examination of the multijet m ff variable 
with varying jet requirements should reveal the general mass scale for the low-lying strongly coupled superpartners, 
and this information should assist in optimizing jet-p^ requirements on signals with lower statistics. 

In Fig. Q we investigate the effect of changing these jet-pr requirements for Chargino 1 (left panel) and Stau 6 
(right panel). Here we include the signatures denoted multi-jets200 and multi-jetslOO in Table (II]) as well as our 
softer signature multi-jets40 from Table §V\ . In this case we show the signal superimposed on the SM background 
separately to distinguish among the three cases. Loosening the jet requirements populates the lower energy bins in 
m e ff for both the signal and the background, leaving the higher energy bins mostly unaffected. Therefore these softer 
cuts tend to boost the signal significance of models such as Chargino 1 with softer jets, while actually reducing the 
significance of models like Stau 6 where jets carry much more px on average. This is also reflected in the signal 
significances at 1 fb _1 in Table §V\ . For these models all three signature definitions yield a discovery, but for other 
models such as Gluino 2 the softer jet requirement is crucial to making an early discovery. In Fig. ( 10 ) we plot the 
generalized Ht 400 variable for the three models for which this is the leading discovery channel: Chargino 5, Stau 6, 
and Gluino 2. For the last case, this channel produces a discovery in approximately 300 pb _1 while other channels 
must wait for more data for a five-sigma excess. Despite the very light gluino in model Gluino 2, the compressed 
spectrum of this model produces extremely soft jets and a relatively low jet multiplicity for gluino pair production 
events. Expanding the definition to include all visible decay products produces a stronger signal. 

As mentioned above, the peak in the jet effective mass distribution - or the energy at which the signal begins 
to exceed the SM background - is a relatively good indicator of the rough mass scale of the lightest 5£/(3)-charged 
superpartner |49j . Apart from this rough estimate, one might ask if there are other properties of the superpartner 
spectrum that can be "measured," even at such low statistics. We believe that the answer is in principle yes, 
depending on the model. The most well-known measurement technique is to form the invariant mass of opposite-sign 
dilcpton pairs and look for an edge, or end point, in the distribution. This procedure is often performed after a 
flavor subtraction is done (i.e., the combination e + e~ + — e + /i~ — e~/i + is formed) to reduce SM and SUSY 

combinatorial backgrounds. After 1 fb _1 there is unlikely to be sufficient numbers of OS dilepton events to perform 
this subtraction, or even to identify a true edge. Even after 2 fb _1 the statistics are sufficiently low to make this 
measurement difficult in all but the most favorable model points. We illustrate this in Fig. (Ill, normalized to 2 fb _1 
of data. It is tempting to see an edge in the invariant mass distribution for Chargino 1 (see inset for a clear view). 
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Generalized Ht Signature 
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FIG. 10: Distribution of Ht, defined by Eq. Q but applied to the hardest four visible objects in the event as opposed to the 
four hardest jets with an isolated lepton veto. The signal plus background distribution is given for the three benchmarks for 
which this was the best discovery channel: Chargino 5 heavy (red) line, Stau 6 thin (green) line and Gluino 2 medium (blue) 
line. The component of the distribution made up of SM background events is indicated by the hatched region. 
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FIG. 11: The opposite-sign dilepton invariant mass is displayed for events satisfying the signature definitions denoted OS- 
dileptons40 in Table pT|. The models shown are Chargino 1 heavy (red) line, Chargino 4 thin (green) line and Stau 3 medium 
(blue) line. As with previous distributions the curves are for the signal + SM background, with the background component 
indicated by the hatched region. The inset shows the distribution for Chargino 1 alone, to exhibit the edge at approximately 
50 GeV. 



But the low statistics are evidenced by the scale on the vertical axis. The sharpness of the edge is largely an artifact of 
simulating our signal sets at the 5 fb _1 level and then rescaling to a smaller integrated luminosity. Nevertheless, this 
edge is indeed real and accurately describes the mass distribution between the lightest and second-lighte st n eutralino. 
The first hint of the "spoiler mode" ~~ ^ Z Xi present in the Stau 3 benchmark is also present in Fig. (11). 

Potentially more promising are measurements based on the invariant mass distribution of events with precisely two 
b-tagged jets. In Fig. ( 12 ) we give di-bjet invariant mass distributions for four models with particularly strong signals 
in this channel. In the right panel we give two models for which the b-jet channels are the strongest discovery modes, 
Gluino 3 and Stopl, and with a sufficient amount of data it should be able to measure the mass difference between 
the gluino and the lightest neutralino. 
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Invariant Mass Distribution of B-Jets 
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FIG. 12: The invariant mass distribution of b-jet pairs for events with precisely two b-tagged jets, each with > 40 GeV, is 
given for events which satisfy our initial cuts of transverse sphericity St > 0.2 and at least 200 GeV of f: T . No lepton veto is 
imposed on these distributions. As with previous distributions the curves are for signal + SM background, with the background 
component indicated by the hatched region. Left panel: Chargino 1 heavy (red) line and Chargino 5 dash-dotted (blue) line. 
Right panel: Gluino 3 heavy (red) line and Stop 1 dash-dotted (blue) line. 



VI. CONCLUSIONS 



We have given here an analysis of the potential of the LHC in early runs at y/s — 7 TeV to discover supersymmetry 
with 1 fb _1 of data. We have carried out an independent analysis of the standard model backgrounds at y/s = 7 
TeV and 1 fb _1 of integrated luminosity which are generally consistent with a previous study |4]. Our analysis is 
done within the framework of the MSSM. However, the parameter space for soft breaking for the MSSM is rather 
large (with over a hundred parameters) and thus intractable. This parameter space is significantly reduced in the 
models we consider. Thus in the analysis of this work we use the framework of mSUGRA which has four parameters 
and the sign of /i, as well as supergravity models with nonuniversality in the gaugino sector which increases the 
number of parameters to six and the sign of \x. The analysis is done under the imposition of theoretical constraints 
which include the radiative breaking of the electroweak symmetry, R parity conservation, and conservation of charge 
and color, as well as under the experimental constraints which include constraints on the relic density from WMAP, 
constraints from flavor changing neutral currents, and the experimental lower limit constraints on the masses of 
the Higgs bosons and on the sparticle masses. The residual model space which passes both the theoretical and the 
experimental constraints as described above is still large and we use a selection criteria in its further exploration. 
Noting that the central theme of the analysis is the discovery of supersymmetry in early runs, our search for models 
is then narrowed to such models as are discoverable with y/s = 7 TeV and 1-2 fb^ 1 of data. Within this general 
theme we choose the models as broadly as possible to encompass as many diverse possibilities as possible. Specifically 
we select models so that all allowed NLSPs such as the chargino, the stau, the gluino, the CP odd Higgs, and the 
stop are included. Using the above criteria we have presented a set of benchmarks for early discovery in Table (III). 
We observe that in Table (III) all the NLSPs mentioned above are represented. Further, in Table (III) the range of 
inputs varies widely. Thus, for example, mo ranges from 101 to 2225, m^M ranges from 313 to 755, Aq ranges from 
—2531 to 2710 (all masses in GeV), and tan/3 ranges from 5.7 to 47.2. The analysis of this work shows that an mi 
mass up to 400 GeV and mo up to 2 TeV could be accessible for certain combinations of soft parameters. Further, 
the precise set of discovery modes in which the benchmarks will become visible are identified. It is shown that most 
of the benchmarks have at least a minimum of five discovery channels and several of them have as many as 10 in 
which the SUSY signal is discoverable. We have also exhibited the dependence of the number of discovery channels as 
a function of the integrated luminosity which shows its rapid increase as the integrated luminosity increases toward 
1 fb^ 1 for the set of discoverable models investigated. Thus if SUSY is found in one of the signature channels for 
any of the benchmarks given in Table (III), then it should also show up in other signature channels as identified 
in Fig. ([5]) providing important cross-checks for the discovery of supersymmetry. All the benchmarks exhibited are 
consistent with the current limits on the spin independent neutralino-proton cross section from XENON100 and 
CDMS-II dark matter detectors. Further, all of the benchmarks will be accessible in the next generation of dark 
matter experiments. 
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VII. TABLES 



A display of signatures/cuts used in early discovery analysis at the LHC at -^s = 7 TeV 
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TABLE II: List of signatures and cuts used in the early discovery analysis. Our notation is as follows: £ = e, /x, n(x) is the 
number of object x in the event, and PT(x n ) is the transverse momentum of the n th hardest object x. For the case of Pt(t) 
we take this to mean the visible part of the pr from a hadronically decaying tau. We required $ T > 200 GeV and a minimum 
transverse sphericity of 0.2. The symbol | should be read as the logic "or": i.e. the cut n (r + \ t~) = 2 would be read "the 
number of r + equals 2 or the number of t~ equals 2." 
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Benchmarks for Early Discovery at -^/s = 7 TeV with 2 fb 
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TABLE III: Benchmarks for models discoverable at the LHC at y 7 ^ = 7 TeV with 2 fb -1 of integrated luminosity. The model 
inputs are given at Mgut = 2 x 10 16 GeV, sign(^i) = +1, and Si = 0. The displayed masses are in GeV. All models satisfy 
REWSB and the experimental constraints as discussed in Sec. ffify . The spin independent direct detection cross section, asi, is 
exhibited as well as the cross section ctsusy for the production of supersymmetric particles at yfs = 7 TeV. Our analysis shows 
that all the models listed in this table are discoverable at the 5a level above the background in several channels as exhibited 
in Fig. (5}. 



Sparticle mass spectra for the benchmarks 
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f\ 


fx 


fx 


ti 


ti 


o"susy (pb) 


0.3 


1.4 


0.4 


0.5 


0.6 


0.4 


2.0 


0.2 


2.3 


3.8 


M 


630 


425 


426 


301 


393 


443 


388 


432 


1213 


691 


m N 1 


297 


157 


181 


188 


241 


280 


159 


211 


232 


324 


m N 2 


630 


198 


408 


303 


395 


415 


189 


408 


508 


689 


m cf 
m fl 


629 


198 


405 


295 


390 


409 


187 


406 


508 


687 


1456 


167 


194 


192 


248 


289 


176 


221 


1546 


1436 


m h 


1032 


541 


504 


529 


329 


506 


497 


615 


258 


357 


m-g 


594 


771 


835 


817 


834 


817 


709 


913 


532 


422 



TABLE IV: An exhibition of the light sparticles for the benchmarks given. These benchmarks are listed in Table ( |III[ ). All 
masses are given in GeV. 
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Signature Name 




SM 


CI 


C4 


C5 


G2 


G3 


S3 


S6 


Tl 


Multi-jets200 


Events 


47 


91 


68 


105 


28 


16 


49 


88 


12 


S/-/B 


13.3 


9.9 


15.2 


4.1 


2.4 


7.2 


12.7 


1.8 


Multi-jetslOO 


Events 


180 


401 


225 


213 


114 


83 


77 


171 


69 


s/Vb 






















29.9 


16.8 


15.9 


8.5 


6.2 


5.8 


12.8 


5.2 


Multi-jets40 


Events 


215 


497 


316 


218 


135 


107 


77 


176 


85 


S/VB 




33.9 


21.6 


14.9 


9.2 


7.3 


5.3 


12.0 


5.8 


#t400 


Events 


965 


1035 


501 


496 


286 


183 


156 


419 


143 


s/Vb 




33.3 


16.1 


16.0 


9.2 


5.9 


5.0 


13.5 


4.6 


Multi-bjetsl 


Events 


188 


460 


188 


175 


51 


126 


50 


102 


86 


s/Vb 




33.5 


13.7 


12.8 


3.7 


9.2 


3.6 


7.5 


6.3 


Multi-bjets2 


Events 


46 


157 


49 


69 


7 


57 


19 


39 


45 


s/Vb 






23.1 


7.3 


10.1 




8.4 


2.8 


5.7 


6.6 


l-lepton40 


Events 


367 


45 


20 


74 








30 


38 


27 


s/Vb 






















2.4 


1.0 


3.9 






1.6 


2.0 


1.4 



TABLE V: LHC discovery channels after 1 fb -1 of integrated luminosity for selected benchmark models from Table (IVl. All 
signatures require transverse sphericity St > 0.2 and at least 200 GeV of $ T . For each signature the number of signal events 
is given, as well as the signal significance if there are sufficient signal events. We have also included a much weaker multijet 
signature (multijets40) in which the four jets are all required merely to satisfy p^, > 40 GeV. This signature also appears in 
Fig. Q. 
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